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Abstract 

We calculate the model-independent radiative corrections to the Dalitz plot of Kfi decays to the order 
of (a/ir)(q/Mi), where q is the momentum transfer and M% is the mass of the kaon. The final results 
are presented, first, with the triple integration over the variables of the bremsstrahlung photon ready to be 
performed numerically and, second, in an analytical form. These two forms are useful to cross-check on 
one another and with other calculations. This paper is organized to make it accessible and reliable in the 
analysis of the Dalitz plot of precision experiments and is not compromised to fixing the form factors at 
predetermined values. It is assumed that the real photons are kinematically discriminated. Otherwise, our 
results have a general model-independent applicability. 

PACS numbers: 14.40.Df, 13.20.Eb, 13.40.Ks 
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I. INTRODUCTION 



Analyses of the low-energy s — > u and d — >■ u semileptonic transitions play a decisive role 
in our understanding of the interplay between weak and strong interactions and the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix. Current determinations of \V u d\ and |\4 S | 
provide the most precise constraints on the size of the CKM matrix elements and yield the most 
stringent test of CKM unitarity. In particular, the best determination of | V us | is achieved from kaon 
semileptonic (Ki 3 ) and leptonic (Kn) decays, whereas the most precise determination of \V u d\ is 
obtained from superallowed + — >■ + Fermi transitions and also, to a minor extent, from baryon 
semileptonic decays -most commonly from neutron beta decay- and pion beta decay. 

Over the past years, precise measurements have been made in both kaon and baryon 
semileptonic decays [1]. In these experiments the statistical errors are rather small, and more 
effort has been put into the reduction of the systematic errors, which are mainly of two kinds. The 
first one comes from the different shortcomings of the experimental devices. The second one, of a 
theoretical nature, comprises assumptions about form factors and radiative corrections. 

On the one hand, while the leptonic part of s — > u and d — >■ u semileptonic transitions is 
unambiguous, the hadronic part is affected by flavor SU(3) symmetry breaking in the form factors. 
For Ki 3 decays there is a remarkable simplification because only the vector part of the weak current 
has a nonvanishing contribution and only two form factors appear. These form factors are protected 
by the Ademollo-Gatto theorem [2] against SU(3) breaking corrections to lowest order in (m s —rh) 
so that the theoretical approach to compute them is under reasonable control within the limits of 
experimental precision. In contrast, baryon semileptonic decays play an ancillary role due to the 
presence of both vector and axial- vector currents and the appearance of six form factors. Although 
the leading vector form factor is also protected by the Ademollo-Gatto theorem, large theoretical 
uncertainties arise from first-order SU(3) breaking effects in the axial-vector form factors. 

On the other hand, the precise and reliable calculation of radiative corrections to the various 
measurable quantities accessible to current experiments is a rather difficult task. Despite the 
enormous progress achieved in the understanding of the fundamental interactions with the standard 
model, no first-principles calculation of radiative corrections to kaon and baryon semileptonic 
decays is yet possible. Thus, in dealing with these radiative corrections it is very important to have 
an organized systematic approach to deal with the complications and technical difficulties that 
accompany them. These radiative corrections depend on the details of strong and weak interactions 
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present in the decay vertex; i.e., they have a model-dependent part. They also depend on the charge 
assignments of the decaying and emitted hadrons. Furthermore, their final form depends on the 
observed kinematical and angular variables and on certain experimental conditions. 

Various articles on radiative corrections to K l3 decays exist in the literature, each with a 
different emphasis. Some of the earliest attempts can be traced back to the works b y G insberg 
1 3-^], Becherrawy |7], Garcia and Maya [8], and more recently by Cirigliano et.al. [9, 10], Bytev 
et. al. Jllll . and Andre [12], to name but a few. Ginsberg, for instance, focused on the calculation 
of radiative corrections to the lepton spectrum, Dalitz plot, and decay rates of Kiz decays, by 
assuming a phenomenological weak K-n vertex; the results depend on a cutoff. Becherrawy, in 
turn, used a particular model of the strong interactions, whereas Garcia and Maya extended to M/ 3 
the procedure proposed by Sirlin ft 1 311 . originally introduced to study the radiative corrections to 
the charged lepton spectrum in neutron beta decay. Cirigliano et. al. computed these radiative 
corrections in the framework of the chiral perturbation theory and accounted for virtual photons 
and leptons. Bytev et. al., in turn, presented an alternative approach to remove the ultraviolet 
cutoff dependence and set it equal to the W mass. Finally, Andre focused his analysis on Kf 3 
decays and included contributions from outside the kinematically allowed three-body region of 
the Dalitz plot. 

As for baryon semileptonic decays, following the analysis of Sirlin of the virtual radiative 
corrections in neutron beta decay ft 1 311 and armed with the theorem of Low for the bremsstrahlung 
radiative corrections [14, Jj|, it has been shown that the model dependence that arises in this latter 
case contributes to the order of (a/n)(q/MB) n for n > 2 and that such model dependence can be 
absorbed into the already existing form factors [16], while the model-independent contributions 
can be extracted from orders n = 0,1. Here M B is the mass of the decaying baryon, and q 
is the four-momentum transfer. This approach has been implemented to compute high-precision 
radiative corrections to various observables in baryon semileptonic decays (decay rates and angular 
spin- asymmetry coefficients) by considering unpolarized and polarized decaying and emitted 
baryons [| 1 VL 1 1 8fl - for all possible charge assignments. 

In this paper, we reexamine the calculation of radiative corrections to decays up to the order 
of (a/7r)(g/M 1 ), where Mi denotes the mass of the kaon. Our analysis builds on earlier works, 
particularly Ref. [8], but unlike it, here we will focus on the radiative corrections to the Dalitz plot 
of Ifjg decays. Our approach to the calculation of these radiative corrections has been to advance 
results which can be established as much as possible once and for all. This task is considerably 



3 



biased by the experimental precision attained in given experiments and by the available phase 
space in each decay. The radiative corrections obtained to the order of approximation considered 
here are then suitable for model-independent experimental analyses and are valid to an acceptable 
degree of precision. We assume that the form factors can be extracted from experiment, and thus 
we do not consider them here. Their model-independent parts contain information only on QED. 

Our results will be presented in two forms: one where the triple integration over the real photon 
three-momentum is left indicated and ready to be performed numerically and another one, of 
analytical nature, where such an integration has been performed. Both forms can be used to 
numerically cross-check one another. However, the analytical result, although tedious to feed into 
a Monte Carlo program, can help in the reduction of computer time because the triple integration 
involved does not have to be performed within the Monte Carlo calculation every time the energies 
of the charged lepton and the pion or the form factors are changed. 

This paper in organized as follows. In Sec. HU we introduce our notation and conventions; 
also we briefly discuss some basic aspects of the phenomenology of kaon semileptonic decays 
and define the boundaries of the allowed kinematical region, usually referred to as the Dalitz 
plot. We will specialize our calculation to the three-body region of this plot. We proceed to 
analyze in Sec. [Ill] the virtual radiative corrections and properly identify the infrared-divergent 
term. In Sec.HVj we discuss the bremsstrahlung radiative corrections and also extract the infrared 
divergence and the finite terms that accompany it. We are left with several triple integrals, mainly 
over the photon variables, which at any rate can be computed numerically. We, however, proceed 
further in order to analytically perform these integrals, and we present the resultant expressions. 
We also present some cross-checks for completeness. We summarize our findings in Sec.|Vj where 
we also carry out a detailed numerical evaluation of these radiative corrections at several points 
of the Dalitz plot; this way we are able to compare our results with others already published, in 
particular, with Refs. |9|,LL0|]. We close this paper with a brief discussion of our results and some 
concluding remarks in Sec.|VTl This work is complemented with two appendixes. In Appendix lAl 
we discuss some features on how to deal with the bremsstrahlung differential decay rate in order to 
properly isolate the infrared divergence. In AppendixE we list the triple integrals over the photon 
variables that emerge in our calculation. 
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II. BASICS ON KAON SEMILEPTONIC DECAYS 



In this section, we provide a basic description of the phenomenology about kaon semileptonic 
decays in order to introduce our notation and conventions. 

For definiteness, we concentrate on the analysis of the semileptonic decay of a positively 
charged kaon, namely, 

K + ( Pl )^7r (p 2 )+£+(l) + p e (pl), (1) 

because the corresponding charge conjugate mode can be analyzed in a similar way. The four- 
momenta and masses of the K + , n°, £ + , and vi will be denoted by p\ = (Ei, pi), p 2 = (E 2 , p 2 ), 
I = (E, 1), and p v = (E®, p°), and by Mi, M 2 , m, and m u , respectively. Because no assumptions 
will be made about the size of m compared to M 1; our results will be valid and applicable to both 
Kf 3 and K^ 3 decays. Without loss of generality, the reference system we will use throughout 
this paper is the center-of-mass frame, which corresponds to the rest frame of K + , so that all 
expressions which are not manifestly covariant will apply to this frame. In this regard, quantities 
like P2, I, or p u will also denote the magnitudes of the corresponding three-momenta, unless 
explicitly stated otherwise. In passing, let us mention that the direction of a generic three-vector p 
will be denoted by a unit vector p. 

The allowed kinematical region in the variables E and E 2 for kaon semileptonic decays, usually 
referred to as the Dalitz plot, is determined by energy and momentum conservation. This region is 
depicted in Fig. Q3 in which we have distinguished two areas labeled as TBR (three-body region) 
and FBR (four-body region). The former (the shaded region in this figure) is bounded by 

E 2 in <E 2 < £™ x , m<E<E m , (2) 

where 

pimax,min 

2 V ^ X " ' 2(M X -E±l) 
and 

E m = ^l-(M 1 2 -M 2 2 + m 2 ), (4) 

whereas the latter is delimited by 

M 2 <E 2 < £™ n , m<E<E c , (5) 

where 

1 77l 2 

E c = -(Mi - M 2 ) + — r. (6) 

c 2 V 1 2J 2(Mi - M 2 ) 



1 M 2 

£max,nnn = L( Ml -E±l) + , (3) 



E 2 




FIG. 1: Kinematical region for Kj^ decays as a function of the energies of the charged lepton and the 
emitted pion, E and E2, respectively, for (a) and (b) decays. In this work, areas TBR and FBR are 
loosely referred to as the three- and four-body regions of the Dalitz plot, respectively. 

The distinction between these two areas is subtle and requires some discussion. First, let us 
notice that in order to find an event with energies E and E 2 in area FBR a fourth particle must exist 
(in our case, such a particle will be a photon) which will carry away finite energy and momentum. 
In contrast, in area TBR this particle may or may not do so. It turns out that area FBR is exclusively 
a four-body region, whereas area TBR is both a three- and a four-body region. Hereafter, areas 
TBR and FBR will be loosely referred to as the three- and four-body regions of the Dalitz plot, 
respectively. 

We can now proceed to the construction of the uncorrected transition amplitude (i.e., the 
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amplitude without radiative corrections) for process (OQ). By neglecting scalar and tensor 
contributions and assuming that only the vector current contributes, such an amplitude, denoted 
here by M , has the structure 

M = C K ^=V: s W,( Pl ,p 2 ) {uMOMl)} > ( 7 ) 

where Gp is the Fermi constant as extracted from muon decay, V us is the relevant CKM matrix 
element, and C# = l/y/2 is a Clebsch-Gordan coefficient. In order to avoid cumbersome 
expressions, the factors Gf, V us and Ck will be included into a single factor G hereafter. 
Furthermore, V£ and u v are the Dirac spinors of the corresponding particles, M = 7 M (1 + 75), 
and the metric and 7-matrix convention we adopt is the standard one (see, for instance, Ref. 1190), 
except that our 75 has the opposite sign. 

The hadronic matrix element W^{p\^p 2 ), on the other hand, can be written as 

W fl {p 1 ,p 2 ) = (n°(p 2 )\u 1 ,s\K + (p 1 )) 

= U(q 2 )(pi+P2), + f-(q 2 )(pi-P2),- (8) 

Here q = p\ — p 2 is the four-momentum transfer and f±(q 2 ) are dimensionless form factors which 
are relatively real if time reversal holds. 

Let us remark that the momentum transfer dependence of the form factors has been discussed 
extensively in the literature; in particular, Ref. [1] provides an overview about the current situation 
of both the theoretical and experimental bent. In summary, a common practice advocated in K^ 3 
analyses is the use of a linear parametrization such as 



f±(g 2 ) = /±(o) 



(9) 



'2 J 

where A+are slope parameters. Most data are well described by this parametrization for 
A_ = [1]. Recent analyses, however, have opted to use the (A + , A ) parametrization instead, 

which is based on the introduction of an alternative set of form factors, namely, 

■2 



W) = hid 2 ) + »^ m H<l% (10) 



M'{-Ml 

where the vector and scalar form factors f+(q 2 ) and fo(q 2 ) represent the p and s wave projections 
of the crossed channel (0\u'J^s\Kti), respectively. Again, if it is assumed that /+ follows a linear 
parametrization and /_ is constant, then 

fo(q 2 ) = /o(0) f 



1 + Ao Mf. 



(11) 
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Finally, another parametrization which has gained a renewal of interest is the dispersive 
parametrization [20], which is based on dispersive techniques and the low-energy K-n phases 
to parametrize the form factors, both scalar and vector. 

As we have pointed out in the introductory section, we assume that the form factors are 
determined from experiment and will not be determined here. Therefore, once the decay amplitude 
M is defined the uncorrected differential decay rate for process (OQ), represented here by dT , can 
readily be obtained by standard techniques. The differential decay rate is thus given by IU9I1 

1 d 3 p2 m d 3 l m v d 3 p° v 



dTr 



2M 1 2£ 2 (2tt) 3 E (2vr) 3 E° (2tt) : 



■(27r) 4 5 4 ( Pl -P2-/-^)^|M c 



(12) 



spins 



The integral over the three-momentum of the neutrino is straightforward. Besides, the integrals 
over the angular variables facilitate if, without loss of generality, we orient the coordinate axes in 
such a way that £ + is emitted along the +z axis and n° is emitted in the first or fourth quadrant of 
the (x, z) plane. Then, the only nontrivial angular integration is over 6> 2 , namely 

1 



dTr 



(2tt) 3 2M 



1 — 1 spins 



where y = cos 6> 2 and 



yo 



E° 

V 



vl 



2p 2 l 



(13) 



(14) 



After some rearrangements, the resultant expression can be cast into the compact form 



dr = Aq dO, 



(15) 



where A is a function of the kinematical variables and depends quadratically on the form factors. 
It can be organized as 



Ao = 4 0) l/ + (? 2 )| 2 + < J Re [/ + (g 2 )/!(g 2 )] + A^f^q 



|(0) 



i Co) | 



(16) 



with 



(0) 



-4 + 



8E, 



4M 9 2 1QE 
2 + 



M 1 Ml M x 



1 - 



E, 



E 



Mi Mi 



+ 



m 



Ml 



6E 2 M 9 2 8E 
M. Mf M. 



Af 



2m 
~Mj 



2 r 



3- 



2E 2 
Mi 



Ml 



2 1 



AE m 
Ml Mi + Ml 



and 



.4 



(0) 



m 



2E, 



+ 



Ml 



m 

2 772 



Ml I Mi Ml Ml 



m 

Ml_ 



7) 



(18) 



(19) 
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and the factor dQ reads 

dQ = c2 G f\Vus\ M 3 dEdE2 _ (20) 

A glance at expressions (fT7T)-(fT9l) reveals why experiments about decays usually determine 
both f + and /_, whereas experiments about K e3 decays are only sensitive to / + because /_ 
comes along with a term proportional to the positron (electron) mass squared, which renders 
its contribution negligible. In Tables U and [XT] we present numerical evaluations of at 
various points (E, E 2 ) of the Dalitz plots of K^ 3 and K^ 3 decays, respectively, for the sake of 
completeness. 

Now that we have some insight into the main features of kaon semileptonic decays, we proceed 
to the calculation of radiative corrections. Let us first discuss the virtual case and later the 
bremsstrahlung one. 



III. VIRTUAL RADIATIVE CORRECTIONS 



To first order in a, the virtual radiative corrections in decay arise from the analysis of 
the Feynman diagrams depicted in Fig. [2] These virtual radiative corrections can be computed 
in complete analogy with the approach implemented by Sirlin II 1 311 in the study of the energy 
spectrum of the electron in neutron beta decay, subsequently extended to the analysis of hyperon 
semileptonic decays by Garcia and Juarez In this approach, the virtual radiative corrections 
can be separated in two parts. One part is model-independent and finite in the ultraviolet region 
and fully contains the infrared divergence. The other one depends on the details of the strong 
interactions. All in all, the separation relies on general principles such as Lorentz covariance, 
analyticity of the strong and weak interactions, and the validity of QED. Garcia and Maya [8] have 
already implemented the procedure to Mj^ decays, so we will partially borrow their methodo 
to achieve our goal. Further details on the procedure can be found in the original papers JsL [l3 , 
so here we limit ourselves to describe only a few salient facts. 

First of all, from the analysis of Fig. [2] one can determine that Figs. I2a)-|2b) comprise the 
graphs in which a photon is emitted from a hadronic line or the intermediate vector boson and 
is absorbed by the charged lepton. To the order of (a/n)(q/Mi), these diagrams yield the 
contribution [|8|] 



G a f ^ u D fia (k)u u O x {-2l a + ^ a )ve 
Mvi ~ 7f« J d k k*-2l-k + ie 



k 2 — 2pi ■ k + ie 

(21) 
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TABLE I: Values of A\ ', Eqs. (fl71)-(fl9i in the TBR of the process K+ -)• vr° + e+ + v e . The entries 
correspond to (a) A[ 0) , (b) x 10 5 , and (c) Af ] x 10 6 . The energies E and E% are given in GeV. 



E 2 \E 


0.0123 


0.0370 


0.0617 


0.0864 


0.1111 


0.1358 


0.1604 


0.1851 


0.2098 


0.2592 


0.0810 


0.3810 


0.6010 


0.7410 


(a) 
0.8010 


0.7810 


0.6810 


0.5010 


0.2410 


0.2468 




0.2110 


0.4510 


0.6110 


0.6910 


0.6910 


0.6110 


0.4510 


0.2110 


0.2345 




0.0410 


0.3010 


0.4810 


0.5810 


0.6010 


0.5410 


0.4010 


0.1810 


0.2222 






0.1510 


0.3510 


0.4710 


0.5110 


0.4710 


0.3510 


0.1510 


0.2098 






0.0010 


0.2210 


0.3610 


0.4210 


0.4010 


0.3010 


0.1210 


0.1975 








0.0910 


0.2510 


0.3310 


0.3310 


0.2510 


0.0910 


0.1851 










0.1410 


0.2410 


0.2610 


0.2010 


0.0610 


0.1728 










0.0310 


0.1510 


0.1910 


0.1510 


0.0310 


0.1604 












0.0610 


0.1210 


0.1010 


0.0010 


0.1481 














0.0510 


0.0510 














(b) 










0.2592 


0.3804 


0.3375 


0.2947 


0.2518 


0.2090 


0.1661 


0.1233 


0.0804 


0.0376 


0.2468 




0.3483 


0.3054 


0.2625 


0.2197 


0.1768 


0.1340 


0.0911 


0.0483 


0.2345 




0.3590 


0.3161 


0.2733 


0.2304 


0.1875 


0.1447 


0.1018 


0.0590 


0.2222 






0.3268 


0.2840 


0.2411 


0.1983 


0.1554 


0.1126 


0.0697 


0.2098 






0.3375 


0.2947 


0.2518 


0.2090 


0.1661 


0.1233 


0.0804 


0.1975 








0.3054 


0.2625 


0.2197 


0.1768 


0.1340 


0.0911 


0.1851 










0.2733 


0.2304 


0.1875 


0.1447 


0.1018 


0.1728 










0.2840 


0.2411 


0.1983 


0.1554 


0.1126 


0.1604 












0.2518 


0.2090 


0.1661 


0.1233 


0.1481 














0.2197 


0.1768 




0.2592 


0.0265 


0.0265 


0.0265 


0.0265 


(c) 
0.0265 


0.0265 


0.0265 


0.0265 


0.0265 


0.2468 




0.0801 


0.0801 


0.0801 


0.0801 


0.0801 


0.0801 


0.0801 


0.0801 


0.2345 




0.1337 


0.1337 


0.1337 


0.1337 


0.1337 


0.1337 


0.1337 


0.1337 


0.2222 






0.1872 


0.1872 


0.1872 


0.1872 


0.1872 


0.1872 


0.1872 


0.2098 






0.2408 


0.2408 


0.2408 


0.2408 


0.2408 


0.2408 


0.2408 


0.1975 








0.2944 


0.2944 


0.2944 


0.2944 


0.2944 


0.2944 


0.1851 










0.3479 


0.3479 


0.3479 


0.3479 


0.3479 


0.1728 










0.4015 


0.4015 


0.4015 


0.4015 


0.4015 


0.1604 












0.4551 


0.4551 


0.4551 


0.4551 


0.1481 














0.5087 


0.5087 





where k and D^ik) are the photon four-momentum and propagator, respectively. The model 
dependence that arises fan these diagrams is contained in Tyfio* explicit form is no, needed 
here but can be obtained without difficulty by following Ref. 1131]. We only point out that is 
regular as k — > and transverse in the sense that ^T Am = 0. 

In the same way, the graph of Fig. (2c) contains the positron wave function renormalization, 
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TABLE II: Values of A\ ', Eqs. (fI7Tl -(fl9i in the TBR of the process K + -> vr° + fi+ + v^. The entries 
correspond to (a) a[ 0) , (b) , and (c) x 10. The energies E and E2 are given in GeV. 



E 2 \E 


0.1131 


0.1280 


0.1429 


0.1578 


0.1727 


0.1876 


0.2025 


0.2174 


0.2322 


0.2480 










(a) 
0.6777 


0.5767 


0.4466 


0.2874 


0.0990 


0.2361 




0.7065 


0.7045 


0.6734 


0.6131 


0.5238 


0.4053 


0.2577 


0.0810 


0.2242 


0.5682 


0.6070 


0.6166 


0.5972 


0.5486 


0.4708 


0.3640 


0.2280 


0.0629 


0.2123 


0.4571 


0.5075 


0.5288 


0.5210 


0.4840 


0.4179 


0.3227 


0.1984 


0.0449 


0.2004 


0.3460 


0.4080 


0.4410 


0.4447 


0.4194 


0.3650 


0.2814 


0.1687 


0.0269 


0.1885 




0.3086 


0.3531 


0.3685 


0.3549 


0.3120 


0.2401 


0.1390 




0.1766 






0.2653 


0.2923 


0.2903 


0.2591 


0.1988 


0.1094 




0.1647 






0.1774 


0.2161 


0.2257 


0.2062 


0.1575 


0.0797 




0.1528 








0.1399 


0.1612 


0.1532 


0.1162 


0.0501 




0.1409 












0.1003 


0.0749 
















(b) 










0.2480 










0.0520 


0.0409 


0.0298 


0.0188 


0.0077 


0.2361 




0.0895 


0.0785 


0.0674 


0.0564 


0.0453 


0.0343 


0.0232 


0.0122 


0.2242 


0.1050 


0.0940 


0.0829 


0.0718 


0.0608 


0.0497 


0.0387 


0.0276 


0.0166 


0.2123 


0.1094 


0.0984 


0.0873 


0.0763 


0.0652 


0.0542 


0.0431 


0.0320 


0.0210 


0.2004 


0.1138 


0.1028 


0.0917 


0.0807 


0.0696 


0.0586 


0.0475 


0.0365 


0.0254 


0.1885 




0.1072 


0.0961 


0.0851 


0.0740 


0.0630 


0.0519 


0.0409 




0.1766 






0.1006 


0.0895 


0.0785 


0.0674 


0.0563 


0.0453 




0.1647 






0.1050 


0.0939 


0.0829 


0.0718 


0.0608 


0.0497 




0.1528 








0.0983 


0.0873 


0.0762 


0.0652 


0.0541 




0.1409 












0.0807 


0.0696 






0.2480 










(c) 
0.0110 


0.0110 


0.0110 


0.0110 


0.0110 


0.2361 




0.0331 


0.0331 


0.0331 


0.0331 


0.0331 


0.0331 


0.0331 


0.0331 


0.2242 


0.0552 


0.0552 


0.0552 


0.0552 


0.0552 


0.0552 


0.0552 


0.0552 


0.0552 


0.2123 


0.0773 


0.0773 


0.0773 


0.0773 


0.0773 


0.0773 


0.0773 


0.0773 


0.0773 


0.2004 


0.0994 


0.0994 


0.0994 


0.0994 


0.0994 


0.0994 


0.0994 


0.0994 


0.0994 


0.1885 




0.1215 


0.1215 


0.1215 


0.1215 


0.1215 


0.1215 


0.1215 




0.1766 






0.1435 


0.1435 


0.1435 


0.1435 


0.1435 


0.1435 




0.1647 






0.1656 


0.1656 


0.1656 


0.1656 


0.1656 


0.1656 




0.1528 








0.1877 


0.1877 


0.1877 


0.1877 


0.1877 




0.1409 












0.2098 


0.2098 







and, after mass renormalization, it contributes [8] 

G a 



M 



v 2 



W x (pi,p 2 ) / tfkD^ikfaO 



4tn fv^n (-/ + ™)(2fo + 7aW/(2f M +M (22) 



A" 



2m 2 (k 2 + 21 ■ k + is) 2 



a/2 87T 3 2 

which is infrared-divergent. 

Finally, Figs. Wid) and[2te) contain the graphs in which the photon is emitted by a hadronic line 
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(e) 

FIG. 2: Feynman diagrams, to first order in a, which yield virtual radiative corrections in decay. 
The wavy, broken, and continuous lines represent virtual photons, pseudoscalar mesons, and fermions, 
respectively. The blobs represent the effects of the strong interactions, and, at the weak vertex, they also 
represent the effects of details of the weak interactions. 
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or the intermediate vector boson and is absorbed by the same hadronic line or another one or the 
intermediate boson. These figures thus yield [80 



M 



v 3 



G a w / \- n /"j4/n /nf^i-^Pi-^a . M/ 
W\(p 1 ,p 2 )u v O\Ve I a kD^k)—^ — ^ - | + M V3 



M C y 3 + M'y 3 , 



(k 2 — 2pi-k + is) 2 



(23) 



where M v , explicitly defined in Eq. (1231) . is infrared-divergent whereas M' v , which can be written 



'v 3 

in the most general way as 



[oi(g 2 )(pi +P2)n + fl2(g 2 )(pi -P2) m ]^7m( 1 +1b)vi, 



(24) 



1/3 tt'^ 1 

is infrared-convergent. The explicit expressions of the additional form factors a^(g 2 ) introduced 



in Eq. (1241) are not needed for our purposes. 

At this stage it is already possible to achieve the separation referred to above. The transition 
amplitude with virtual radiative corrections to the order of (a/7r)(q/Mi) is constituted by two 
parts My and My. The former is model-independent and gauge-invariant and contains in full the 
infrared divergence, whereas the latter contains all the model dependence. M y is given by adding 
together the term proportional to the first summand within square brackets in Eq. (|2~TT) . My 2 , and 
My 3 . After some analysis we find 



MV = M 07 ^$ 1 (£) + M Pl ^-$ 2 (£), 

where the functions $>i(E) and & 2 {E) are given by 
1 



(25) 







arctanh (3 



In 


"A" 






m 



— — (arctanh (5) 2 + —L 



h - b 2 




hi - b 2 


[ -b 2 _ 


-v 


l-b 2 _ 



+ — arctanh (3 


~M 2 - 


M 1 E(l + f3 2 )' 


+ 


"3 m 2 ' 


In 


"Mi" 








m 


-im 


ri-bi 

l-b 2 _ 


jln 


~M{ 

m 


— arctanh (3^ 
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>_ T' 







and 



where 



1-/3 2 



/3 



—arctanh (3 



1 + 



M\E - m 2 



Mi 



m 



h,2 



m 



M\E ± M x l 



r2 i _2 



M{ + m 2 - 2M 1 E. 



(26) 
(27) 

(28a) 
(28b) 
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Here (3 = l/E, L is the Spence function, and A is the infrared-divergent cutoff. This divergent 
term will be canceled by its counterpart in the bremsstrahlung contribution to be discussed in the 
next section. 

Similarly, the second term in Eq. (|25l) is given by 

M " = -i&S w ^° xM - <29) 

As for the model-dependent part My, it is given by the sum of the term proportional to 
the second summand within square brackets in Eq. (f2~Tj) and M'y 3 . The sheer impossibility of 
computing analytically the integrals over the photon four-momentum involved in My leads us to 
implement Lorentz invariance instead, so we get [8Q 

cx G 

My = 7=[ai(q 2 ,p+-l)(pi+P2) t t + a2(q 2 } p + -l)(p 1 - p 2 )^u v ^(l + 75)^, (30) 

where a" and a!^ are some other form factors which now depend on q 2 and p + ■ I = (pi + p 2 ) ■ I. 
A close inspection of My reveals that it has the same structure as M [Eq. ©], so one is prompted 
to absorb the former into the latter by redefining the original form factors f±, namely, 

M' = M + My 

= ^=[f + (q 2 ,p + ■ Z)(pi +P2)„ + fUq\P+ ■ l)(pi-P2)»]u v (p v )O ll Vz{l), (31) 

where the modified form factors f, and f'_ have a new dependence in the positron and emitted 
pion energies other than the ones in the q 2 dependence of the original form factors. By rearranging 
some terms of the order of 0(a), these form factors can be written generically as 

f + (q\p+ ■ I) = U(q 2 ) + -a+( P+ ■ I), (32a) 

7T 

f'_(q 2 , P+ ■ I) = f-(q 2 ) + -a^(p + ■ I), (32b) 

7T 

where a + and a_ are some other functions which contain all the model dependence and the prime 
on M in Eq. (1311) will be used as a reminder of this fact. 

Let us remark that the introduction of the modified form factors f±(q 2 , p+ ■ I) in Eq. (1311) is not 
a withdrawal of the approach. Ultimately, these modified form factors are the ones which can be 
experimentally measured and are the ones that provide information about the strong interactions. 

By gathering together partial results, the transition amplitude for decay with virtual 
radiative corrections is then given by 



My = M 



a: 



2tt 
14 



$!(£)] +M P1 ^$ 2 (£), (33) 



where some terms of second order in a have been rearranged in the above equation to stress the 
fact that to first order in a only the modified form factors f±(q 2 ,p+ • I) appear in it. 

Armed with the transition amplitude with virtual radiative corrections, My, we can obtain the 
corresponding differential decay rate dT v by means of a long and tedious but otherwise standard 
procedure. Assuming for definiteness complex form factors, the final expression can be cast into 



dT v = dQ \A' (l + + A'y 

V 7T / V7T 



where 



(34) 



A' v = A<P\f4(?,p + -l)\ 2 + 4 y) Re [f + {q\p + • 0£V,P+ ' 01 + A^\fUq 2 ,P + • 0l a , (35) 



with 



A iv) - — 

Al ~ My 



AE 



M| 



A 



(V) 



Ml 

AE 
~M~i 



E 
Mi 



1 + 



M| 



2 Mi 
2E 2 



M| 



+ 



m 



Ml Mi 



2M iy / 2Mf 

~ mT 



E2 

Mi 



~Mf" 



and 



(V) 



£ 2 
2Mf 



2E 2 M| m 2 (Mi - E 2 ) 
1 ~ M7 + M 2 MfE 



(36) 
(37) 

(38) 



Equation (|34l ) is our first partial result for decays. It is the differential decay rate with 
virtual radiative corrections to the order of (a/n)(q/Mi). It contains an infrared-divergent term 
in &\(E) which at any rate will be canceled when the bremsstrahlung radiative corrections are 
added. Let us stress the fact that all the model dependence arising from the virtual contribution 
has been absorbed into the modified form factors (|32l) which are the ones that enter into Eq. (l34l) . 
The primes in this equation are an indicator of this fact. Let us now discuss the bremsstrahlung 
contribution. 



IV. BREMSSTRAHLUNG RADIATIVE CORRECTIONS 

In this section we now turn to the analysis of the emission of a real photon in the process 

K+( Pl ) 7r°(p 2 ) + £ + (l) + u e (p v ) + 7 (A0, (39) 

where K + denotes a positively charged kaon and n° a neutral pion, whereas £ stands for a 
positively charged lepton (£ = e + or and v t its accompanying neutrino. We again point 
out that the charge conjugate mode of process (l39l) can be analyzed likewise. Here, 7 represents 
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a real photon with four-momentum k = (u, k). As in the case of virtual radiative corrections, we 
will use the rest system of the decaying kaon. Accordingly, energy and momentum conservation 
yield M% = E + E 2 + E v + to and = p 2 + 1 + + k, where the neutrino energy and momentum 
in the presence of the photon are, respectively, 



EL = E°, 



Vv = Vv 



(40) 



where E„ and pj; are the energy and three-momentum of the neutrino in the nonradiative process 
(OQ), respectively. 

A quick glance at Fig. Q] reveals that the TBR of the Dalitz plot is the region where the three- 
body decay CD and the four-body decay (1391 ) overlap completely. In contrast, the FBR is where 
in process (|39~1) neither of the energies of the neutrino and the real photon can be zero. Therefore, 
strictly speaking, process (1391) is kinematically allowed to occur anywhere in the joined area TBRU 
FBR of Fig. [TJ If we assume that real photons can be discriminated in an experimental setup, then 
our analysis of bremsstrahlung radiative corrections will consider process (l39l ) restricted to the 
TBR. 

To first order in a, the bremsstrahlung radiative corrections in decays arise from the 
analysis of the Feynman diagrams depicted in Fig. [3] We now need to obtain the transition 
amplitude of process (|39l) . M B , keeping all the terms of the order of (a/ it) (q/M-\ ) explicitly. 
This can be achieved in a model-independent way by using the Low theorem 11141. Il5ll . which 
asserts that the radiative amplitudes of orders 1/k and (k)° can be determined in terms of the 
nonradiative amplitude without further structure dependence. Therefore, the direct application of 
the Low theorem allows us to express M B as 

4 

M b = J2 M b v (41) 



»=i 



where 



M 



eG 
72 



Wx[u v O x vA 



l-e pi-e 



l-k 



M 



B 2 



Pi ■ k 

V£, 



M Ba 



— (f 



f- 



Pi-e 
Pi ■ k 



A e A 



u u O\Ve, 



and 



M Bi = ^ 



eG 

7! 



Pi ■ e 
Pi ■ k 



q ■ k — q ■ e 



dW x _ n 
-p^p u v (J\Ve. 



(42) 
(43) 
(44) 

(45) 
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i K 

(a) 



J K 

(b) 




(d) 



FIG. 3: Feynman diagrams which yield bremsstrahlung radiative corrections in decay. The wavy lines 
represent real photons, the broken lines represent pseudoscalar mesons, and the continuous lines represent 
fermions. 

Here e is the polarization four-vector of the photon. Let us notice that, strictly speaking, M b a in 
Eq. (|45l) will contribute to the order of 0(q 2 ) to the decay rate so it will be suppressed with respect 
to the others and can be ignored in our analysis. 

^fow, the differential decay rate with bremsstrahlung radiative corrections dT B can be written 



In order to correctly account for the contribution of the unobserved photons to the Dalitz plot of 
process d39l) restricted to the TBR, we need to perform a careful analysis to delimit the integrations 
over the kinematical variables in (|46l) . We have already mentioned that the system of reference 
we use is the rest frame of the decaying kaon. Accordingly, with all generality, we use the same 




dT B 



1 1 mm, 



d 3 p 2 d 3 ld 3 p u d 3 kS 4 ( Pl -p2-l-pu~k) \ M b\ 2 - (46) 



(2tt) 8 2M x 4:E 2 EE„cj 



spins, e 
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orientation of the coordinate axes as in the previous cases: namely, the direction of emission 
of £ + coincides with the +z axis and 7r° is emitted in the first or fourth quadrant of the (x, z) 
plane. With this choice we are left with five out of 12 variables of the final state. Two of them 
should be, of course, the energies E and E 2 of £ + and it , respectively, whose limits are given in 
Eq. ©. The other three nontrivial variables can be grouped into two sets, namely, (k, cos 9 k , (p k ) 
and (cos #2, cos9 k , Lp k ), where k, 9 k , and ip k are the magnitude of the three-momentum and the 
polar and azimuthal angles of the photon, respectively, and 9 2 is the polar angle of n°. The former 
set is more suitable for dealing with the infrared-divergent terms, whereas the latter is more useful 
in the analysis of the infrared-convergent contributions. For this latter case one has p 2 l = cos 9 2 = 
y, 1 • k = cos 9 k = x, and p 2 • k = cos 9 2 cos 9 k + sin 9 2 sin 9 k cos (p k , and the photon energy can 
be expressed as 



where 



w = 2lr (47) 



F = 2p 2 l(y -y), (48a) 
D = El + Ix + p 2 ■ k, (48b) 

and the scalar product p 2 • k can readily be written in terms of x and y. 
As for dT B we consider it convenient to organize it as 

dv B = dv Bl + dv B2 + dr B3 + dv Bi + dr B5 , (49) 

where dT Bl ~ E|M Bl | 2 , dT B2 ~ ^ [|M Ba | 2 + 2Re[M Bl M^ 2 ]] , dT B3 ~ ^2Re[M Bl M^ 3 ], 
dY BA ~E2Re[M B2 My,and dT Bs ~ J2\M Bs \ 2 . 

The contribution dT Bl deserves particular attention because it not only contains an infrared- 
divergent term but also finite ones that come along with it, which must be properly identified. For 
this purpose, let us notice that | M Bl \ 2 can be conveniently separated as 



l2 e 2 G 2 M? , „ , 
t-^ 2 mm v ^-^ 



; n 2 

I ■ k p\- k 



(50) 



spins, e 

where I, pi, k, and e within the square brackets in Eq. (l50l) are understood to be four- vectors. 
Besides, A is given in Eq. (fT6l) . and the function B reads 

B = B^Uiq 2 )] 2 + B 2 Re {f + ( q 2 )fl( q 2 )} + B 3 \f_(q 2 )\ 2 , (51) 



18 



where 



+ 



Mf 



2E 2 

TT + M 2 



+ 



M x 4 Ml 

2Eu(l - fix) 
Ml 

2u(M 1 + E 2 -p 2 -k) 
Ml 



8E 2m 2 
M\ + ~Mf 



Ml 



2E 2 

1 H + 

+ Mi Ml 



i _ ^ + Mi 

Mi 



2w 2u{E 2 - p 2 • k) 
Mi Mf 



4£ m 2 

m7 + Ml 



(52) 



5, 



2F 

Ml 



+ 



4a; 
Mi 



2Sw(l - fix) 
Ml 

2E 2 



rn 
Ml 



2 + Mi 



+ 



2E 2 4E 



-1 + 



M| 



Mi 



2E 

Ml + Mi 



TO 

M 2 



+ 



Ml 

u(E 2 - p 2 ■ k) 
Ml 



+ 



TO 



Mi Mi Mi 2 

"m7 



and 



TO 



+ 



2Eu{\ - fix) 



TO 



M 2 M 2 



(53) 



(54) 



M 2 M 2 M 2 

Equation (|5Q|) has been purposely separated the way it stands because the factor A is precisely 
the one that is needed to cancel the infrared divergence contained in its virtual counterpart Eq. (|34|) . 
Now, the extraction of the infrared divergence and the finite terms that come along with it can be 
conveniently performed by following either one of the approaches implemented by Ginsberg Q3[| 
(introduced in the analysis of Kj^ decays) or by Tun et. al. [21] (implemented in the analysis of 
hyperon semileptonic decays). Ultimately, it has been shown that both approaches are equivalent 
and yield the same results j22\\ . For convenience, in this work we follow the latter approach, and 
in Appendix lAl we briefly describe some important aspects of the procedure. 

On the other hand, we should exercise some caution when performing the pending sum over the 
photon polarization in Eq. (l42l) . which also raises an important issue. In all infrared-convergent 
terms the ordinary covariant summation can be used, namely, J^(e • u)(e ■ v) = —u ■ v, where 
u = (uo, u) and v = (v ,v) are arbitrary four-vectors and tu = k, with k the magnitude of 
k. However, in the infrared-divergent terms the longitudinal degree of polarization of the photon 
must be accounted for. This can be carried out by using the Coester representation Q23[| in which 



v^/ w \ (u-k)(v-k) 
2> ' M )( e ' y ) = u ' v 772 



(55) 



where lu 2 = k 2 + X 2 and A is a fictitious mass given to the photon to regularize the infrared 
divergence. 
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By taking into account all the above arguments, cITb 1 can be written as 

dv Bl = dr Bi + dv ic Bv 



(56) 



where dT B , the piece containing the infrared divergence, can be evaluated through a direct 
application of Eq. (IA9I ); the resultant expression is 



dT 



'2tt 



-dn^p 2 A 



dx 



k 2 x 2 /u 2 



2tt rki £.2 J 

k 2 x 2 /co 2 



l JO 

3n/2 pk 2 fc2 ^ 

I dip k \ dk—\g~(9 2 ^ r 

tt/2 Jk 4 U 

3tt/2 r-k 2 fc2 

- I dip k \ dk—\g + (6 2/l 

tt/2_ Jk 4 {u-/3kx) 



(u) — (3kx) 2 
l-k 2 x 2 /u 2 



where we have used the definition 121(1 



9 ± (0 2 ) 



sin 6^ 



a sin 62 



b cos 62 ' 



(57) 



(58) 



and the factors that appear in Eqs. ((571) and (1581) are defined in Appendix El 

Following the lines of Ref. Kill , we can split the range of integration over k into two intervals, 
namely (0, Ak) and (A/c, fcj), with i = 2,4. This approach is general enough to avoid taking 
the limit Ak — > in the infrared-divergent terms because Ak cancels exactly before performing 
the tp k integration. In contrast, the limit Ak — > can be taken in the infrared-convergent terms 
because they are regular in Ak. Let us notice that in the interval < k < Ak, the function g~(6* 2 ) 
can be expanded in powers of k; the expansion yields j2l\ \ 

1 



9 



+ 0(k). 



(59) 



Thus, after some rearrangements, dT Bi becomes 



7T ZTT 



dx 



2t 



, f Ak „k 2 l- k 2 x 2 /u 2 

d<p k / dk — — r Y 

2p2< Jo io w - pkx) 1 



it/2 



</0 
fc4 fc2 y 

dk — \g~{@2^ 1 



Ak 



id 



dip k 

-7T/2 

/.3tt/2 rk 2 y,2 

+ / c?(y9 fc / dk — \g 

Jtt/2 JAk U 

r3n/2 rk 2 ^2 

+ / d(p k / dk— \g 

Jtt/2 Jk 4 



k x jur 





- (5kx) 2 




k 2 x 2 /lo 2 




- (3kx) 2 




k 2 x 2 /co 2 




- (3kx) 2 



+ 0(Ak) 



On the other hand, y can also be expanded in powers of k 112 ill , namely, 

y = y -fk + O(k 2 ), 

20 



(60) 



(61) 



where 



f' = — [El + (p 2 y + l)x + p 2 (l - 2/^) 1/2 (l - x 2 ) 1 / 2 cos <p k ] . 



(62) 



Thus, the last three summands in Eq. (l60l) can be cast into the form (I A 131) . replacing the upper 
limit over the y integration with y — f'Ak, and also by using the definition of the photon energy, 
Eq. (|47T) . namely, 



dT 1T 



^dn p ^A 

7r 2ir 



+ 



2- 



l&il Jo 
1 

dx 
i (1 



d(f k 



1 



x 



dx / dk— — — 

! J u {lo - pkx) z 

r 2w nyo-f'Ak 1 



(63) 



2^2^ i-i (l - fay j j_ x yo-y. 

Now, the (pk integration in the first summand within the square brackets in Eq. (1631) can be trivially 
performed, whereas in the second summand we have rewritten the integral with the x integration 
outermost so we can easily compute the y integration. This results in 



dT 



a 



Bi 



7T 



dVLAo 



2 

— — I dx 

4:71 



1 , f Ak „ k 2 1 - k 2 x 2 /co 2 
x Jo u (u - pkxy 



X 



(3x) 



2- 



difk In 



f'Ak 
l + 2/o 



(64) 



Finally, the two integrals in Eq. d64l make up the function defined in Eq. (96) of Ref. [21], 
which contains an infrared-divergent piece with the correct coefficient to cancel its counterpart in 
the virtual radiative corrections. With this last result, we have achieved the proper identification 
of the infrared divergence and all the finite terms that come along with it. 

We have already pointed out that the approach introduced in Ref. [4] to deal with the infrared 
divergence is equivalent to the one described here. The equivalence is established through the 
relation Q x — Io(E, E 2 ), where Iq(E, E 2 ), given in Eq. (27) of Ref. [4], in our notation reads 



h{E, E 2 ) 



arctanh/3 



2 In 



m 
X 



2 In 
-In 



+ In 



Vn 



mri 2 m 



4(E + l)rj 



-v 



4r_, 



P 



L 



4r_ 



2mE®(q 2 — m 2 



(65) 



where 

(E + l)r ± 
and 



m 



-t 2 E 



± 



-i 2 



t 2 = 7] m (Ap 2 l - 7) m ), 



EJ {T ~ m 2 ) - -t 2 E 



2p 2 l(l + y ). 



1/2 



- -m 2 t 4 
16 



(66) 



(67) 
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In summary, the final form of dT B reduces to 



a 



dT% = -dQAoOt. (68) 

7T 

As for the infrared-convergent contributions, the use of the set of variables (cos 9 2 , cos 9 k, <fk), 
corresponding to the form (IA13I) . allows us to discern a group of triple integrals which make up 
dT B , dTB 2 , ■ ■ ■ , dTB 5 - Skipping details, these contributions read 

dT Bi = -dSl-^ [(Ai + A 2 + A 3 )|/ + (g 2 )| 2 + (-A 2 - 2A 3 )Re [f + (q 2 )fl(q 2 )] + A 3 |/-(? 2 )| 2 ] , 



71 M( 



(69) 



dY B2 = -dn—t [(A 4 + A 5 + A 6 )|/+(g 2 )| 2 + (-A 5 - 2A 6 )Re [f + (q 2 )fl(q 2 )] + A 6 |/_(g 2 )| 2 ] 



71 M( 



(70) 



dT B3 = -dn-j (A 7 + A 8 )|/ + (g 2 )| 2 + A 7 Re[/ + (g 2 )/:(g 2 )] - A 8 |/_(g 2 )| 2 ] , (71) 

7T M 1 

dv Bi = -dn^ [(A 9 + a 10 )|/ + (? 2 )I 2 + A 9 Re if + (q 2 )f-(q 2 )} - Aiol/-(? 2 )l 2 ] , (72) 

dT Bb = -dQ-^ [A u \f + (q 2 )\ 2 + 2A n Re [f + (q 2 )fl(q 2 )] + An|/-(? 2 )| 2 ] • (73) 

Here Aj (i = 1, . . . , 11) are functions of E and E 2 and contain triple integrals over the relevant 
angular variables. They are listed in Appendix |B] for completeness. 

Equations (1681) and (T69b-(T73T) constitute the first partial result for the bremsstrahlung radiative 
corrections. The infrared divergence and the finite terms that come along with it have been properly 
identified. As a result, the bremsstrahlung differential decay rate can be obtained by performing 
numerically the triple integrals that make up the functions Aj at various points (E, E 2 ) of the Dalitz 
plot. This could require a considerable effort in an actual analysis. 

We can proceed further, however, and also compute analytically the integrals displayed in 
Eqs. (I69l)-(l73l). The main aim of having a fully analytical expression available relies on the 
fact that it can be useful in the reduction of computing time in a Monte Carlo simulation of 
an experimental setup. For this task, then, a considerable amount of work can be saved if we 
realize that such integrals can be expressed in terms of the 9i functions originally introduced in the 



and 



analysis of bremsstrahlung radiative corrections in hyperon semileptonic decays 1121 , 
For conciseness, one has 



22, 



24, 



25] 



Ai = ^ (EY 3 - 2EY 2 + 29 ) , (74) 



^ 2 -2M[ 



-29 + iz x + (E + El)Y 2 - E°Y 3 
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(75) 



A, 



P2IE 



[{2p 2 ly - m 2 )Y 3 - 2p 2 l^Y l - 2Z 2 + 2E(1 - $ 



(76) 



A, 



P2P 



(E° u - 2E)[9 7 - 2E{9 3 - 4 )] + (3EE° U - p 2 ly )9 3 - 3E 



73 — (74 



4 - 



+ 2E Vo - E(l - /3 2 ) [2£°0 2 - 6 + 2£(0 2 - 3 )] - EE% - £0 O - ^9 



(77) 
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2Afj 



771 

[m 2 (E + £j) - E(l 2 + EEl)] 9 3 - — 7 + E(l 2 + p 2 Zy - m 2 )9 A 

El 

+ E(E- E° u )r] + E 2 9 + Pp 2 m 2 9 u - El 2 9 10 + —9 14 + Em 2 Y 2 - E( n 



(78) 



Ae = M 2 m2 (P^yo-l 2 -EEy 3 + E(m 2 E° u ~2Ep 2 ly )9 A + EE J(2W 10 -9 
+ ^[2EE° U - 2l 2 +p 2 l{l - y )]r]o - Em 2 9 + 2(E 2 ( W - m 2 ( n ) 



(79) 



A, 



P2I 
2Mi 



[£(1 - /3 2 )(£ + £°) + + / 2 ] 3 - {m 2 + 2££°)0 4 - E(9 - Vo ) 



E 



-(2l/39 10 - /30 14 ) - -{I - P 2 )9 7 + C11 - Z^0 6 



(80) 



and 



A 



A* 



P2I 
AM 2 



(y (1 - 2/o) - /) fyo + (^2/2/0 + £> 2 )#4 - ^Cio + m 2 9 



1 



+ 2Z 2 p 2 (y o 05 - Y x ) - m 2 [El + /3(l + p 2 y )} 9 3 + -E°J(2W 1Q - 9 



14 J 



A 9 = TTTT (P2^0^4 - - Clo) 



2M X 



A 



10 



V2IE 
AM 2 



Kvo - P2ly o 4 + |£v(2Z0i O - 0u) + C10 



11 



Ep 2 l 
8M 2 



f3E° u (2W 10 - U ) + \ [2E° U - 2/3/ + 0p 2 (l - y Q )] Vo + 2(3p 2 l(y 9 5 - Y 1 ) 



(81) 

(82) 
(83) 

(84) 



where the secondary functions 7] , Q m , Z r , and Y s , given in Ref. [25], depend on the 0, functions. 

Finally, a further simplification of Eqs. (T7"4l)-(l84"l) yields a fairly compact form of the 
bremsstrahlung differential decay rate. It can be organized as 



a 



dT B = -dn (A 9 1 + A B ) 

7T 



(85) 
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where 



.4 



B 



A B) \U 



2\ (2 



+ A { 2 B) Re[f + ( q 2 )fl(q 2 )] + A^ , \f- 



( R) 

and the functions AS read 



M l Am 



2E 2 5E 



m 
' E 2 

2 

m 



2EU 
3 



M{ 
i ~ 

m 
2E 2 



F 2 

mF (1 + ^ 



+ 



m 

.U, ' W 
2E 9 E AE 2 



2M 2K 
El - 3E - p 2 f3y 
3E(E + E° v )+p 2 ly 



yo) 



+ 



+ 



4 



M 2 



Mi Mi Mi 4Mi 



3E 2 AE 

9 H H 

+ Mi Mi 



. ,2E 2p 2 ly 
h Mi Ml 
APE 2 P2 l 



M 2 



Mi 



m 

AM 2 
2p 2 l 
' Mi 
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ifc + ^6 + 



-2 + 







,2\|2 



2E 

M 



(86) 



3 m 2 
AM 2 



7 7 



013 + Mf ^ 14 



£ MiE 

T y l9 + 772*20, 



2£ 



M 2 



Mi 



(87) 



Mf 1 
m 2 4^2^ 



.4 



1 + ™ 
'2V + ^ 2 



3 m 
2 + £Mi 



Eo 
+ — 



2Mi 



3 + 



2E + E 2 
2Mi 



+ 



2Mi 



-0 



5- 



and 



.4 



p 2 lm 2 
Mf 



l-9 + (Mi - £ 2 + /3p 2 yo)0 3 - (Mi - £ 2 )#4 - ^ 5 



#4 

(88) 



(89) 



Expression (f8~5l) constitutes our second partial result for the bremsstrahlung differential decay 
rate of Kj^ decays, restricted to the TBR of the Dalitz plot. In spite of the lengthy expressions it 
consists of, its final form is easy to handle. Both infrared-divergent and infrared-convergent terms 
that appear in it have been rigorously identified and appropriately extracted. Let us now discuss 
the series of cross-checks performed on the functions A*. 



A. Numerical cross-checks 

At this stage, we consider pertinent to cross-check the functions Aj in their analytical forms 
(T7"4l)-(l84l) against their counterparts (IB1|) - (IB1 1|) with the triple integrals indicated, in order to 
ensure that our results are correct. These cross-checks consist in computing numerically the 
triple integrals and then contrasting the outputs with the direct evaluations of the corresponding 
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TABLE III: Values of J (see the text) in the TBR of the process K + — > tt° + e + + v e by (a) performing 
numerically the triple integration in Eq. d90T > and (b) evaluating straightforwardly Eq. d91~| ). The energies E 



and E 2 are given in GeV and J is dimensionless. 



E 2 \E 


0.0123 


0.0370 


0.0617 


0.0864 


0.1111 


0.1358 


0.1604 


0.1851 


0.2098 












(a) 










0.2592 


0.3793 


-0.5527 


-2.0511 


-3.6155 


-4.9209 


-5.7042 


-5.7307 


-4.7722 


-2.5782 


0.2468 




0.7234 


-0.8941 


-2.6557 


-4.1663 


-5.1285 


-5.2902 


-4.4211 


-2.2974 


0.2345 




2.0511 


0.3903 


-1.5012 


-3.1741 


-4.3068 


-4.6362 


-3.9305 


-1.9761 


0.2222 






1.6926 


-0.3135 


-2.1392 


-3.4412 


-3.9464 


-3.4197 


-1.6509 


0.2098 






2.9937 


0.8784 


-1.0965 


-2.5672 


-3.2501 


-2.9063 


-1.3252 


0.1975 








2.0634 


-0.0593 


-1.6975 


-2.5578 


-2.3960 


-0.9997 


0.1851 










0.9642 


-0.8405 


-1.8758 


-1.8919 


-0.6739 


0.1728 










1.9645 


-0.0048 


-1.2106 


-1.3972 


-0.3461 


0.1604 












0.7955 


-0.5726 


-0.9158 


-0.0113 


0.1481 














0.0153 


-0.4537 














(b) 










0.2592 


0.3793 


-0.5527 


-2.0511 


-3.6155 


-4.9209 


-5.7042 


-5.7307 


-4.7720 


-2.5780 


0.2468 




0.7234 


-0.8941 


-2.6557 


-4.1663 


-5.1285 


-5.2901 


-4.4210 


-2.2971 


0.2345 




2.0511 


0.3903 


-1.5012 


-3.1741 


-4.3068 


-4.6362 


-3.9304 


-1.9760 


0.2222 






1.6926 


-0.3135 


-2.1392 


-3.4412 


-3.9463 


-3.4196 


-1.6507 


0.2098 






2.9937 


0.8784 


-1.0965 


-2.5672 


-3.2501 


-2.9062 


-1.3251 


0.1975 








2.0634 


-0.0593 


-1.6975 


-2.5578 


-2.3959 


-0.9996 


0.1851 










0.9642 


-0.8405 


-1.8757 


-1.8919 


-0.6739 


0.1728 










1.9645 


-0.0048 


-1.2106 


-1.3972 


-0.3461 


0.1604 












0.7954 


-0.5726 


-0.9158 


-0.0113 


0.1481 














0.0153 


-0.4536 





analytical expressions at various points (E, E 2 ) of the Dalitz plot. As an example, let us define the 
quantity 



P2 i a2 f 1 A i - x 2 f 2 * A fy° F B 1 



where B\ is given by Eq. (1521) . Its analytical counterpart is 

J=^(Ai + A2 + A 3 ), (91) 

with A r above given by Eqs. (|74l)-(|76l). We present in Table UlTl the evaluation of expressions 
( |90| ) and (|9TT) for Kf z decays. We observe a very good agreement, entry by entry, between the 
two forms. Of course, we have also cross-checked the other Aj functions and found very good 
agreements, too. We have repeated the exercise for K^ 3 decays and also found an excellent match, 
so there is no need to reproduce these outputs here. 
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V. FINAL RESULTS AND COMPARISON WITH OTHER CALCULATIONS 



The differential decay rate of Kj^ decays in the variables E and E 2 , restricted to the TBR of 
the Dalitz plot and including radiative corrections to the order of (a/7r)(g/M 1 ), is given by 

dF(K±) = dT v + dT B . (92) 

dT v is given by Eq. (|34l) . For dT B two forms are available. In the first one the triple integration 
over the real photon variables remains to be performed numerically. It is given by the sum of 
Eqs. (|68l) to (|73l) . which are written in terms of the functions Aj, Eqs. (|B 1|) - (1B 1 1|) . listed in 
Appendix [Bj We should point out, however, that the infrared divergence cancels exactly in the 
sum in Eq. (|92| ). The second form for dTs is given by Eq. (f85l) and is expressed in terms of the 
fully analytical functions Aj, Eqs. (I74l)-(l84l). Let us remark, however, that obtaining this latter 



form was feasible by 



semileptonic decays 1121 , 



using resu 



22, 



24, 



ts introduced in the analysis of radiative corrections in hyperon 



25] 



Our analytical result can thus be cast into the compact form 

G% I K/.r 1 2 



dT(K±) = C 2 K ^^MldEdE 2 [A' + ^A' 



(93) 

where A' comprises all the various contributions arising from radiative corrections, namely, 



A' = A'^ + O^ + A'y^ + A'z 

= A 1 \f,(q 2 ,p + ■ 0| 2 + A 2 Re [f + {q\ P+ ■ l)fL*(q 2 ,p + ■ I)) + A 3 \f__(q 2 ,P + ■ (94) 

where the functions Aj are implicitly defined in Eq. (|94l) . Here A' Q , A' v , and A' B , given by Eqs. (fT6l) . 
(1351) . and (f86l) . respectively, are functions of E and E 2 and, to a very good approximation, depend 
quadratically on the modified form factors (l32l) . and the primes on them are an indicator of this 
fact. Besides, $i and <J> 2 are given by Eqs. (|26l) and (1271) . For 9i we can use either Eq. (96) of 
Ref. |21j or Eq. (1651) of the present work since they are equivalent. 

We are now in a position of producing some numerical evaluations which, at the same time, 
will allow us to compare our outputs with the ones obtained within other approaches. In the 
introductory remarks, we pointed out that some treatments about radiative corrections in 



decays are available in the literature J3412D. In the present paper, we put emphasis on the radiative 



corrections to the Dalitz plot of Kj^ decays so a suitable comparison can readily be performed 



with the results presented in Refs. la 1100 for this particular case 
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InRef. H the spin- averaged decay distribution p(y, z) for K13 decays is analyzed, where, in 
our notation, the variables y and z correspond to y = 2E/M\ and z = 2E 2 /Mi, respectively. 
Also, the uncorrected distribution p^\y^z), Eq. (3.5) of Ref. [90, is expressed in terms of the 
kinematical densities Af\y,z), A 2 (y^z), and A^\y,z), Eq. (3.8) of this reference, which in 
turn are written in terms of the variables r t = m 2 /M 2 and r n = M| /M 2 . 

The uncorrected distribution p(°>(y,z) corresponds to our dT (E,E 2 ) of Eq. (fT5l) . and the 
amplitudes Af\y, z) are algebraically equivalent to our A ( f\E, E 2 ) given in Eqs. (fT7T)-([T9l). A 
numerical evaluation of A^ for decay is presented in Table 3 of Ref. Iioll . which can be 
compared with Table Ufa) of the present paper. The agreement, entry by entry, is evident. 1 

On the other hand, the analysis of radiative corrections to Ki 3 decays presented in Refs. ipl 
[Toll was performed in the context of the chiral perturbation theory to the order of 0(p 4 , (m^ — 
m u )p 2 , e 2 p 2 ), with the inclusion of the photon and light leptons as active degrees of freedom. In 
summary, the density distribution p(y, z) with radiative corrections, which is equivalent to dT of 
Eq. ([93]), is written as Q9Q 



p(y,z)=AfS EW (M p ,M z ) [A 1 \f + (t,v)\ 2 + A 2 [f + (t,v)f_(q 2 )]+A 3 \f_(t,v)\ 2 }, (95) 

where the effective form factors f± depend on t = (pi — p 2 ) 2 and v = (pi — I) 2 , 

\f r»2 G 2 F \V us \ 2 Mf 

M = ° K 128.3 ' (%) 

and Sew is tne short distance enhancement factor. 

The effects of virtual photons are contained in the long distance component T c {v , m 2 , M 2 ; M 7 ) 
of loop amplitudes (which produces infrared and Coulomb singularities); it depends on v and the 
masses of the charged lepton m and the charged meson M and has a logarithmic dependence 
arising from the infrared regulator M 7 . On the other hand, the contributions of real soft photons 
are obtained by virtue of the theorem of Low. Therefore, the kinematical densities to the order of 
a read 

My,z) =A?\y,z) [l + A m (y, z)} + A lB (y, z). (97) 

Here A 1R (y,z) contains T c (v, m 2 , M 2 ; M 7 ) and |A^^„ 1} | 2 , where M 1 = Ml^ + Mj 0) is the 
radiative amplitude and Ail ^ and AiJ^ contain terms of orders 1/k and (A;) , respectively. 



1 We anticipated performing this comparison so we adopted the same notation for the amplitudes and evaluated 



dr (E, E 2 ) at the same values of E and E 2 as in Ref. jlCI. 
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Similarly, A m (y, z) comprises the additional terms of_LM 7 | 2 . From (|971) , the radiative corrections 



10] are 



to the Dalitz plot of K i3 decays obtained in Refs. [|9J, 

My, *) ~ A (0) (l/, z) = A?\y, z)A m (y, z) + A lB (y, z). 
These latter contributions should be equivalent to our (a/n)Ai of Eq. (|94b , which read 

A^^ + ^+A^ + A^' 



(98) 





a ' 


7T 


71 . 


91 and 


(1991) 



(99) 



scope of the present paper. We limit ourselves to performing a numerical comparison of the 
available pieces of information instead. For this purpose, we display a few samples of numerical 
values of the radiative corrections (a/ir)Ai of Eq. (|99l) for both Kf z and K^ 3 decays in Tables 
|IV]and|Vl respectively. As anticipated, for Kf 3 decays, (a/V)/^ and (a/7r)A 3 are negligible, 
but t hey are not so for K^ 3 decays. Additionally, the numerical results displayed in Table 4 of 
Ref. [10] for Ai(y, z) — A x (y, z) of K^ 3 are reproduced in Table |VT1 of the present work to carry 
on a comparison with Table ITvTa). An inspection of these tables shows an overall good agreement 
at the first significant digit over most of the entries and even the agreement at the second digit is 
more evident at the lowest entries of each column. There are small differences, however, and they 
may be explained as due to the different approximations used. 

We can go further and provide a preliminary expression for the decay rate of Kf 3 decays by 
following the lines of Ref. 
factors performed in this reference. 

From Eq. (1931) and in analogy with Eq. (7.5) of Ref. 



T{K t 



9I, making use of the detailed determination of the relevant form 

. I2I we have 

2 /(A H 



r 2 it/ i 2 



128tH 



(100) 



where 



4 

Ml 



dE 



dEo 



(o) 



7T 



h + hi\ + + h 2 \ 2 +J 



(101) 



where the integration limits are given in © and the slope parameter A + = 0.0328 ± 0.0033 has 
been estimated in Ref. ||9j]. Notice that Eq. (| 1 001) does not contain the short distance enhancement 
factor Sew yet. 
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TABLE IV: Radiative correction (a/ir)Ai, Eq. d94), in the TBR of the process K + -)• tt° + e + + f e . The 
entries correspond to (a) (a/ir)Ai x 10, (b) (a/V)^ x 10 6 , and (c) (a/n)As x 10 7 . The energies i? and 
E2 are given in GeV. 



E 2 \E 


0.0123 


0.0370 


0.0617 


0.0864 


0.1111 


0.1358 


0.1604 


0.1851 


0.2098 


0.2592 


0.1533 


0.1880 


0.1462 


0.0668 


(a) 

-0.0286 


-0.1220 


-0.1949 


-0.2246 


-0.1726 


0.2468 




0.1810 


0.1580 


0.0902 


0.0011 


-0.0905 


-0.1654 


-0.1998 


-0.1537 


0.2345 




0.1578 


0.1522 


0.0962 


0.0150 


-0.0718 


-0.1445 


-0.1792 


-0.1354 


0.2222 






0.1429 


0.0989 


0.0261 


-0.0551 


-0.1249 


-0.1590 


-0.1168 


0.2098 






0.1321 


0.1004 


0.0363 


-0.0391 


-0.1056 


-0.1389 


-0.0977 


0.1975 








0.1014 


0.0461 


-0.0233 


-0.0863 


-0.1186 


-0.0779 


0.1851 










0.0558 


-0.0075 


-0.0670 


-0.0979 


-0.0568 


0.1728 










0.0654 


0.0083 


-0.0474 


-0.0769 


-0.0333 


0.1604 












0.0243 


-0.0275 


-0.0550 


-0.0019 


0.1481 














-0.0070 


-0.0316 














(t>) 










0.2592 


0.0540 


0.0227 


-0.0011 


-0.0189 


-0.0315 


-0.0392 


-0.0417 


-0.0383 


-0.0269 


0.2468 




0.0270 


0.0028 


-0.0161 


-0.0301 


-0.0393 


-0.0436 


-0.0422 


-0.0339 


0.2345 




0.0281 


0.0034 


-0.0162 


-0.0311 


-0.0412 


-0.0466 


-0.0468 


-0.0414 


0.2222 






0.0034 


-0.0169 


-0.0324 


-0.0434 


-0.0499 


-0.0516 


-0.0496 


0.2098 






0.0034 


-0.0176 


-0.0339 


-0.0458 


-0.0533 


-0.0567 


-0.0589 


0.1975 








-0.0185 


-0.0356 


-0.0483 


-0.0570 


-0.0623 


-0.0700 


0.1851 










-0.0373 


-0.0510 


-0.0610 


-0.0685 


-0.0844 


0.1728 










-0.0393 


-0.0541 


-0.0655 


-0.0759 


-0.1075 


0.1604 












-0.0576 


-0.0710 


-0.0855 


-0.2223 


0.1481 














-0.0788 


-0.1014 




0.2592 


0.0081 


0.0049 


0.0027 


0.0007 


(c) 

-0.0011 


-0.0031 


-0.0054 


-0.0085 


-0.0143 


0.2468 




0.0189 


0.0122 


0.0064 


0.0007 


-0.0055 


-0.0129 


-0.0232 


-0.0424 


0.2345 




0.0338 


0.0226 


0.0128 


0.0031 


-0.0075 


-0.0202 


-0.0381 


-0.0725 


0.2222 






0.0333 


0.0194 


0.0057 


-0.0094 


-0.0278 


-0.0539 


-0.1056 


0.2098 






0.0440 


0.0260 


0.0081 


-0.0117 


-0.0361 


-0.0709 


-0.1432 


0.1975 








0.0325 


0.0102 


-0.0145 


-0.0452 


-0.0898 


-0.1883 


0.1851 










0.0118 


-0.0181 


-0.0556 


-0.1115 


-0.2473 


0.1728 










0.0126 


-0.0230 


-0.0682 


-0.1377 


-0.3407 


0.1604 












-0.0299 


-0.0844 


-0.1726 


-0.7783 


0.1481 














-0.1090 


-0.2302 





With no radiative corrections, the integral (11011) yields h = 0.0965, h\ = 0.3337, and 
= 0.4618, whereas the inclusion of radiative corrections yields h = 0.0958, hi = 0.3303, 
and h% = 0.4557. In order to compare under the same quotations with Ref. [|9j], we use A + = 0.030, 
so radiative corrections cause a decrease of 0.8% in the decay rate. This has to be compared with 
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TABLE V: Radiative correction (a/ir)Ai, Eq. (|94), in the TBR of the process K + -> ir° + fi + + v^. The 
entries correspond to (a) {a/-n)Ai x 10 2 , (b) (a/V)^ x 10 3 , and (c) (a/vr)^ x 10 3 . The energies E and 
E2 are given in GeV. 



E 2 \E 


0.1131 


0.1280 


0.1429 


0.1578 


0.1727 


0.1876 


0.2025 


0.2174 


0.2322 


0.2480 










(a) 

-0.3524 


-0.2563 


-0.2278 


-0.1927 


-0.1099 


0.2361 




0.0248 


0.0020 


-0.0398 


-0.0880 


-0.1331 


-0.1635 


-0.1626 


-0.0952 


0.2242 


0.0933 


0.0777 


0.0417 


-0.0070 


-0.0607 


-0.1106 


-0.1452 


-0.1473 


-0.0809 


0.2123 


0.0887 


0.0782 


0.0466 


0.0014 


-0.0498 


-0.0980 


-0.1315 


-0.1332 


-0.0661 


0.2004 


0.0695 


0.0670 


0.0423 


0.0029 


-0.0435 


-0.0880 


-0.1190 


-0.1193 


-0.0510 


0.1885 




0.0512 


0.0344 


0.0019 


-0.0389 


-0.0789 


-0.1068 


-0.1053 




0.1766 






0.0248 


-0.0004 


-0.0350 


-0.0702 


-0.0946 


-0.0910 




0.1647 






0.0144 


-0.0032 


-0.0315 


-0.0616 


-0.0824 


-0.0766 




0.1528 








-0.0063 


-0.0280 


-0.0530 


-0.0701 


-0.0651 




0.1409 












-0.0448 


-0.0600 
















(t>) 










0.2480 










-0.2690 


-0.1865 


-0.1568 


-0.1271 


-0.0795 


0.2361 




0.0148 


-0.0315 


-0.0767 


-0.1136 


-0.1398 


-0.1527 


-0.1490 


-0.1263 


0.2242 


0.1431 


0.0611 


-0.0079 


-0.0662 


-0.1136 


-0.1490 


-0.1710 


-0.1779 


-0.1804 


0.2123 


0.1527 


0.0680 


-0.0054 


-0.0686 


-0.1213 


-0.1626 


-0.1915 


-0.2090 


-0.2497 


0.2004 


0.1556 


0.0696 


-0.0068 


-0.0739 


-0.1310 


-0.1775 


-0.2136 


-0.2434 


-0.3671 


0.1885 




0.0701 


-0.0094 


-0.0800 


-0.1415 


-0.1936 


-0.2377 


-0.2830 




0.1766 






-0.0124 


-0.0869 


-0.1531 


-0.2114 


-0.2651 


-0.3326 




0.1647 






-0.0159 


-0.0947 


-0.1662 


-0.2321 


-0.2986 


-0.4054 




0.1528 








-0.1044 


-0.1827 


-0.2590 


-0.3470 


-0.5828 




0.1409 












-0.3075 


-0.4629 






0.2480 










(c) 

-0.0058 


-0.0056 


-0.0064 


-0.0076 


-0.0100 


0.2361 




0.0006 


-0.0017 


-0.0041 


-0.0066 


-0.0093 


-0.0126 


-0.0172 


-0.0282 


0.2242 


0.0087 


0.0048 


0.0011 


-0.0027 


-0.0069 


-0.0116 


-0.0175 


-0.0264 


-0.0495 


0.2123 


0.0144 


0.0097 


0.0047 


-0.0006 


-0.0065 


-0.0134 


-0.0224 


-0.0364 


-0.0773 


0.2004 


0.0208 


0.0151 


0.0088 


0.0019 


-0.0059 


-0.0153 


-0.0277 


-0.0477 


-0.1254 


0.1885 




0.0210 


0.0132 


0.0046 


-0.0054 


-0.0174 


-0.0338 


-0.0611 




0.1766 






0.0176 


0.0071 


-0.0051 


-0.0202 


-0.0412 


-0.0787 




0.1647 






0.0220 


0.0094 


-0.0054 


-0.0241 


-0.0510 


-0.1053 




0.1528 








0.0109 


-0.0070 


-0.0305 


-0.0665 


-0.1706 




0.1409 












-0.0445 


-0.1054 







h = 0.09533, hi = 0.3287, and h 2 = 0.4535 evaluated in this reference, which induces a decrease 
of 1.27%. 
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TABLE VI: Radiative correction [Ai(y,z) - A\ '(y,z)] x 10 presented in Table 4 of Ref. 110] and 
reproduced here for comparison with Table llVf a) (see the text). The energies E and E<i are given in GeV. 



Fn \ F 
JZ/2 


0.0123 


0.0370 


0.0617 


U.UOUI 


0.1111 


0.1358 


0.1604 


0.1851 


0.2098 




1494 


u. luy / 


n 1 1 7 a 


U.UOJ-O 


— u.uu / u 


— u. ijyo 


n 997^ 

— U.ZZ t 




—0 1841 


0.2468 




0.1708 


0.1364 


0.0610 


-0.0320 


-0.1236 


-0.1946 


-0.2213 


-0.1638 


0.2345 




0.1558 


0.1378 


0.0732 


-0.0128 


-0.1006 


-0.1704 


-0.1983 


-0.1440 


0.2222 






0.1356 


0.0821 


0.0036 


-0.0796 


-0.1474 


-0.1758 


-0.1240 


0.2098 






0.1321 


0.0898 


0.0190 


-0.0593 


-0.1248 


-0.1533 


-0.1035 


0.1975 








0.0971 


0.0341 


-0.0392 


-0.1021 


-0.1305 


-0.0822 


0.1851 










0.0490 


-0.0191 


-0.0794 


-0.1075 


-0.0597 


0.1728 










0.0639 


0.0010 


-0.0566 


-0.0841 


-0.0348 


0.1604 












0.0214 


-0.0333 


-0.0598 


-0.0020 


0.1481 














-0.0094 


-0.0340 





VI. CONCLUSIONS 

In this work we have obtained the radiative corrections to the Dalitz plot of Kj^ decays to 
the order of (a/ir)(q/Mi), where q is the momentum transfer and Mi denotes the mass of the 
kaon. We have obtained a fully analytical expression which comprises contributions of both virtual 
and real photons, restricted to the three-body part of the allowed kinematical region. Despite its 
length, the analytical form obtained, Eq. (l93l . is quite simple and organized in a way that is easy 
to deal with. Among other properties, it contains all the terms of the order of (a/ir)(q/Mi), it 
has no infrared divergences, it does not contain an ultraviolet cutoff, and it is not compromised by 
any model dependence of radiative corrections. As argued in Sec. Unl the model dependence is 
absorbed into the already existing form factors by adding a function of p + ■ I to f+(q 2 ) and another 
to f_(q 2 ). This, needless to say, is a theoretical problem and as such should be dealt with like that. 
The usefulness of Eq. (1931) . however, could be better appreciated when incorporated into a Monte 
Carlo simulation, since it may reduce the computational time required by the triple integrals. 

We should emphasize that Eq. (1931) is very useful for processes where the momentum transfer 
is not small and thus cannot be neglected. Thus, it is valid for any decay, whether M be 
n ± , K ± , D ± or even B ± . To first order in q it yields terms of the order of (a/ir)(q/Mi) in 
the radiative corrections. The expected error by the omission of higher order terms is around 
(a/ir)(q/Mi) 2 ps 0.0012 in K ± and D ± decays. Being conservative, if the accompanying factors 
amount to 1 order of magnitude increase, then we can estimate an upper bound to the theoretical 
uncertainty of 1.2%. This should be acceptable with an experimental precision of 2%-3%. We 
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envisage further improvements to our calculation by incorporating the effects of the FBR of the 
Dalitz plot and also the decay of a neutral kaon, which requires an extra effort Ii26ll . 
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Appendix A: The bremsstrahlung differential decay rate and the infrared divergence 



In this appendix, we provide an outline of the method used to extract and iso 



ate the infrared 



divergence in the present analysis. The method has been implemented in Ref. 112 111 to the analysis 
of baryon semileptonic decays so we have borrowed and adapted it to our analysis. Since most of 
the material presented here can be found in Sees. Ill and IV of this reference, we have used the 
same conventions and notation. 



The differential decay rate of process (1391) can be written as [19] 

1 1 mm v 
(2vr) 8 2Mi AE 2 EE v u 



dT B = ——— mmu d 3 p 2 d 3 l d 3 Pu d 3 k 5 4 ( Pl -p2-l-Pu-k) V |M B | 2 . (Al) 



spins, t 

Proper integration over the variables involved in Eq. (|A1I) requires that one chooses an 
orientation of the coordinate axes where the integrand acquires its simplest form. For this purpose, 
we orient the coordinate axes in such a way that 1, the three-momentum of the charged lepton, lies 
along the +z direction, and p 2 , the three-momentum of the neutral pion, lies in the first or fourth 
quadrant of the (x, z) plane. 

Integrating (|A11) over the neutrino three-momentum is an easy matter, so we are left with 



2 

2 



dT B = -r—^-^-—^-dEdE 2 d(- cos 9 2 )d(p 2 dQ t — dkd(- cos 9 k )d(p k 5(f (cos 8 2 )) ^ \M B 

^ v spins. e 

(A2) 

where 9 2 [9 k ] and cp 2 [Lp k ] are the polar and azimuthal angles of the pion [photon], respectively, 

/(cos 9 2 ) = El - u - v/(£° - u) 2 - C + a cos 9 2 + b sin 2 , (A3) 
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and 



C = (E® — tu) 2 — pi — I 2 — k 2 — 2kl cosOk, a = 2p 2 (l + k cos 9 k ), b = p2ksm9 k cos</?fc, 

(A4) 

with El = M x -E- E 2 . 

To proceed further we need to solve the equation / (cos 9 2 ) = to find the zeros of the function 
/(cos 9 2 ) in order to perform the integration over <5(/(cos 9 2 ))d(— cos 9 2 ) in Eq. (|A2I) . The use of 
energy and momentum conservation yields 



± aC ± bVa 2 + b 2 -C 2 
2 a 2 + b 2 



cos Of = = — tt. , (A5a) 



bC T Qv / Q 2 + b 2_Q2 
sin0f = ^— . (A5b) 

Since < 2 < 7r, the signs of cos and sin ^ 2 are fixed unambiguously. 



Armed with Eqs. (|A5I) together with energy and momentum conservation relations, the physical 
values of the photon three-momentum can readily be obtained. First, let k 12 denote the zeros of 
the radicand in Eqs. (IA5I) . and, next, let k 3ti denote the values of k that satisfy C = ±a. The 
physical contents behind the plus and minus signs in the latter condition means that 1 and p 2 are 
parallel and antiparallel, respectively, so the photon and neutrino three-momenta must rearrange 
accordingly to satisfy momentum conservation. We thus find 



and 



where 



p 2 ci - ai&i ± y/{p 2 ci - ai&i) 2 + (b\ - c\)(dl - a\) 
h > 2 = 2{^aJ) ' (A6) 



b 1 -2p 2 (p 2 ±l) 

k 3A = 777 T TTT^ ( A7 ) 

2(a 1 ±p 2 cos9 k ) 



a 1 = E° u + lcos6 k , &! = Ef+p 2 2 -l 2 , Cl = 2p 2 E° u , d\ = p\{\ - sin 2 e fc sin 2 <p k ). 

(A8) 

Throughout a careful and detailed analysis, one finds that k 2 and k± are the only physical values 
of k, whereas k\ and k 5 are unphysical ones. In Fig. 51 we have plotted k 2 and fc 4 as functions 
of 9 k and ip k for decay, at E = 111.1 MeV and £ 2 = 222.2 MeV for definiteness. In this 
figure, the uppermost surface corresponds to k 2 so A; 4 lies just right below it. Although /c 2 ,4 are 
physical values, the condition sin8 2 G [0, 1] strongly constrains the accessible regions to them. 
In other words, in the region where < k < k 4 and < cp k < 2ir, only sin^ and cos^^ are 
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HO k , <p k ) 




FIG. 4: Plot of k 2 and fc 4 as functions of 6 k and yjfc for decay, at -E = 111.1 MeV and E 2 = 222.2 
MeV. The uppermost surface represents k 2 so k^ lies right below it. 

allowed, whereas in the region where k 4 < k < k 2 , both cos Of and sin 6*^ are allowed, but this 
time ifk G ( 7r /2, 3n/2). In Fig. [5l we have plotted k 2 and &4 as functions of (fk for Ok = 7r/3, at 
the same values of (£", i? 2 ) as above. The shaded area depicts the accessible values of k which 
lead to allowed values of cos 2 . 

All the above findings can be gathered together in order to express the bremsstrahlung 
differential decay rate as 



dT, 



1 p 2 l mm v 



dEdE? 



(2tt) 6 2 Mi 

-2n rk 4 

dipk I dk 



+ 



dx 



i 

U3 



asin^ 



+ 



r3n/2 


dip k 


rk2 


k 2 
dk — 


In/2 




lk 4 


OJ 


p3n/2 


d(fk 


rk2 


k 2 
dk — 


In/2 




lk 4 


OJ 




(A9) 



spins, e 



where x = cos 6^ and y = cos0 2 . This form of dT B , in terms of the set of variables (k, 2 , Lfk) 
discussed in Sec.|IVj is suitable for handling the infrared divergence although the involved factors 
in the integrand might complicate the evaluation of the integral to a great extent. There is another 
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2 



In 



FIG. 5: Plot of k% and k^ as functions of cp^ for 9^ = 7r/3, at the same values of [E, E 2 ) as in Fig. [4] 
The continuous and broken lines represent k 2 and k^, respectively. The shaded area displays the accessible 
values of k which lead to allowed values of cos 6 2 . 



set, however, which allows one to handle the infrared-convergent contributions. It is expressed in 
the variables (9 2 , k ,(p k ) through 

F 

(A10) 



*=4 

2D' 



where 



F = bi-2p 2 {p2 + I cos 8 k ), 

D = ai + p 2 (cos9 2 cos 9 k + sm 62 sin. 9k cos 6k) 



(All) 
(A12) 



and after some algebraic manipulation we finally obtain 

1 pj mm v r> r ~" 



dT, 



(2tt) 6 2 M x 



dEdE 2 I dx I d(p k / dy-^ ^ |M 



F 
4D' 



B\ 



spins, e 



(A13) 



This latter form is the one actually used in the analysis of infrared-convergent contributions here. 



Appendix B: Integrals over the photon variables 



The various contributions dY B . are constituted by the integrals listed below. 



P2I 
4vr' 



2/0 r2ir 

dy / dip k 
-1 Jo 



1 - 



E + lx 
D 



(Bl) 
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1-3 



A? = ^r--r-r I dx- — 
4tt M 1 J_, (1 

P2I P 2 



— x 2 r yo r 2n 



E 2 + E°Jx + p 2 ly 



D 



E 



dx- 



1-x 2 



**4M»j_r(i-wj-i dy Jo dLPk 



ya 



2tt 



E(l - 0x)(F - m 2 ) | m2 



(B2) 
(B3) 



P2P f 1 

47T 

+ (2E„ - D) 



dx T^Px~ J dy J d ^ k n{ EE ° ~ l2 -Vily - oo(E + lx) - E V E{\ - px) 
E{l-f3 2 ^ 



u + E - 



1-px 



(B4) 



P2P 1 



f 1 1 f yo f 2w 1 f 

- / dx ~, 5- / d V / d Vk -h ^ " ™ 2 £° + w(p 2 /y + -E 2 + - -E-D) 

1J-1 1 - pa: 7_i Jo -D I 



4tt Mi 

+ E(l - px) [cuD + EE U + l(p 2 y + I + cox)} - 



i + e_ E(i-p 2 y 

OJ Uj(l — Px) 



x [EE v uu(2 - px) + p 2 y + ojx) - EuD] j , 



(B5) 



P2P 1 



4tt AM 2 



l dX T^Px~ J l dy J d ^p[ m2 ( EE » + l2 +P2 l yo) + EF u(l-(3x) 



m 



m 2 Eu(l - Px) + —F - E(l - Px 



(E + to) [F + 2(EE„ + p 2 ly + I 2 + lux)] 



- E v [m 2 + 2Eu(\ - Px)] 



+ 



\ + E_E{l-P 2 ) 



LU u(l — Px 

x [2uE(l - Px){EE v + I 2 +p 2 ly + lux) - m 2 uD] j, 



(B6) 



A 7 = 



V A1 

4tt M 



/l ryo /-27T y ( 

dx dy d Vk - { uE(l + Px) -uD- EE° U + I 2 + p 2 ly 
1 J-i Jo Ef { 

+ Y^pi \- EK[l - ^ - ? x ^ E + u)D - E ° Jx - 12 - P2lyo ^ }' 



(B7) 



P2I 1 



yo 



2tt 



LU 



D 



1 - 



1-px 



[m 2 + 2Eu(l - Px)] 



Ei ;] f I + Z 2 + p 2 Zy - uE{l - pz)] + ^m 2 + EE° U (1 - Px) 

LU{L — px) LU 



E(M 1 - E 2 ) - p 2 ly + (E + lo)D\, 



A fl 



Al 

4tt Mi 



1 J-l Jo 



dtp k u 



1 - 



2£„' 
D 



(B8) 



(B9) 
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A 



P2I 1 



10 



4tt M- 



/1 ryo /"27T 
da: / dy d<p k — [-ED + ^(2 - (3x) + p 2 ly + I 2 + lux] ,(B10) 
1 J-i Jo m 



A 



11 



p 2 l E 



1 pyo 
dx 

1 



dy 

1 Jo 



dy k ^ [(1 - Px)(E° u - u - D) + D] . (Bll) 
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